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AMENDMENTS TO THE CLAIMS 

The listing of the claims will replace all prior versions, and listings, of claims in the 
application: 

LISTING OF CLAIMS : 

1. (currently amended): A method for analyzing multivariate images, comprising: 

a) providing a data matrix D containing measured spectral data, 

b) transforming tho data matr i x D, using a wavo l ot transform, ajjpjyjng. 
at least one wavelet transform to the data matrix D to obtain a transformed data 
matrix D comprising an approximation matrix D^ and a detail matrix D^, 

o) threshold i ng tho wavo l ot co o ff i c i onts of th o transformod data matrix 

4) c) performing an image analysis on the transformod data matr i x t T-te 
nhtn i n a r . pnt i nl l y compressed concentration matr i x - €--aod-a approximation 
matrix D 3 to obtain an approximation concentration matrix C a and a spectral 

shapes matrix S, and 

e) d] computing a concentration matrix C from the spat i a ll y comprossod 

concontrat i on matr i x - 6- approximation concentration matrix C„ , the detail matrix 
D rf . and the spectral shapes matrix S . 

2. (canceled) 

3. (currently amended): The method of Claim 1, wherein the at least one wavelet 
transform comprises a Haar transform. 

4. (canceled) 

5. (canceled) 

6. (currently amended): The method of Claim 1, wherein the image analysis of 
step d) c] comprises an alternating least squares analysis and the spatial l y 

compressed concontrat i on matr i x -e- approximation concentration matrix C a and 
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the spectral shapes matrix S are obtained from a constrained least squares 

7. (currently amended): The method of Claim 6, wherein the alternating least 
squares analysis comprises a transformed non-negativity constraint C a > 0 . 

8. (currently amended): The method of Claim 1 , wherein the computing step e) d) 
comprises: 

computing a detail concentration matrix C d from the detail matrix D rf and 

the spectral shapes matrix S: 

combining the approximation concentration matrix C., and the detail 

concentration matrix C rf to provide a transformed concentration matrix C : and 

applying an inverse wavelet transform to the spat i al l y compr e ss e d 
conc o ntration matr i x ~ 6- transformed concentration matrix C to provide the 
concentration matrix C. 

9. (canceled) 
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10. (previously presented): A method for analyzing multivariate images, 
comprising: 

a) providing a data factor matrix A and a data factor matrix B obtained 
from a factorization of measured spectral data D, 

b) transforming the data factor matrix A, using a wavelet transform, to 

obtain a transformed data factor matrix A , 

c) thresholding the wavelet coefficients of the transformed data factor 

matrix A, 

d) performing an image analysis on the transformed data factor matrix 
A and data factor matrix B to obtain a spatially compressed concentration matrix 
C and a spectral shapes matrix S, and 

e) computing a concentration matrix C from the spatially compressed 
concentration matrix C . 

11. (previously presented): The method of Claim 10, wherein the data factor 
matrix A comprises a total of j blocks of data factors A, and the data factor matrix 
B comprises k blocks of data factors B„ thereby providing a concentration block 
C, in step e), and wherein steps a) through e) are repeated sequentially until the 
concentration matrix C is accumulated blockwise, according to 

c = fc, C 2 - 0 M Cj]. 

12. (original): The method of Claim 10, wherein the wavelet transform comprises 
a Haar transform. 

13. (canceled) 

14. (previously presented): The method of Claim 10, wherein the thresholding 
comprises decimating the detail coefficients. 

15. (previously presented): The method of Claim 10, wherein the image analysis 
of step d) comprises an alternating least squares analysis and the spatially 
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compressed concentration matrix C and the spectral shapes matrix S are 



obtained from a constrained least squares solution of min 

c,s 



AB T -CS T . 



IIF 



16. (original): The method of Claim 15, wherein the alternating least squares 
analysis comprises a transformed non-negativity constraint. 

17. (previously presented): The method of Claim 10, wherein the computing step 
e) comprises applying an inverse wavelet transform to the spatially compressed 

concentration matrix C to provide the concentration matrix C. 

18. (previously presented): The method of Claim 10, wherein the computing step 
e) comprises projecting the product of the data factor matrix A and the data 
factor matrix B from step a) onto the spectral shapes matrix S from step d), 

according to mini AB T - CS T || and subject to appropriate constraints. 



19. (original): The method of Claim 10, wherein the data factor matrix A 
comprises a scores matrix T and the data factor matrix B comprises a loadings 
matrix P, and wherein T and P are obtained from a principal components 
analysis of the measured spectral data D, according to D = TP T . 

20. (original): The method of Claim 19, wherein T and P represent the significant 
components of the principal components. 

21 . (previously presented): The method of Claim 1 , wherein the data matrix D is 
weighted. 

22. (previously presented): The method of Claim 10, wherein the data factor 
matrix A and the data factor matrix B are weighted. 
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23. (new): The method of Claim 1 , wherein the wavelet transform applying step 
b) comprises: 

folding the data matrix D into a (x+1)-dimensional multiway array D 
consisting of x spatial dimensions and 1 spectral. dimension comprising p spectral 
channels, wherein x = 1 , 2, or 3; 

applying an independent wavelet transform to each of the x spatial 
dimensions for each of the p spectral channels to provide a transformed multiway 

array D ; 

partitioning the transformed multiway array D into a multiway 
approximation array D a and a multiway detail array D d ; and 

unfolding the multiway approximation array D a to obtain the 
approximation matrix D a and the detail matrix D d ,. 

24. (new): The method of Claim 23, wherein the computing step d) comprises: 

computing a detail concentration matrix C d from the detail matrix D d and 
the spectral shapes matrix S; 

combining the approximation concentration matrix C a and the detail 

concentration matrix C d to provide a transformed concentration matrix C ; 

folding the transformed concentration matrix C into an (x+1)-dimensional 

transformed concentration matrix C ; 

applying an independent inverse wavelet transform to each of the x spatial 
dimensions of the transformed concentration array C to obtain the multiway 
concentration array C; and 

unfolding the multiway concentration array C to obtain the concentration 
matrix C. 
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25. (new): The method of Claim 1 , wherein D is a 2D data matrix comprising m 
rows and n columns and the wavelet transforms W are applied according to 

(W„®WjxD = D. 

26. (new): A method for analyzing multivariate images, comprising: 

a) providing a data matrix D containing measured spectral data, 

b) applying at least one wavelet transform to the data matrix D to 

obtain a transformed data matrix D comprising an approximation matrix D a , 

c) performing an image analysis on the approximation matrix D a to 

obtain a spectral shapes matrix S, and 

d) . computing a concentration matrix C from the data matrix D, and the 
spectral shapes matrix S. 

27. (new): The method of Claim 26, wherein the at least one wavelet transform 
comprises a Haar transform. 

28. (new): The method of Claim 26, wherein the image analysis of step c) 
comprises an alternating least squares analysis and an approximation 

concentration matrix C a and the spectral shapes matrix S are obtained from a 
constrained least squares solution of min||D a -C a S T 

29. (new): The method of Claim 28, wherein the alternating least squares 
analysis comprises a transformed non-negativity constraint C a > 0 . 

30. (new): The method of Claim 26, wherein the computing step d) comprises 
projecting the data matrix D from step a) onto the spectral shapes matrix S from 

step c), according to min||D-CS T || F , subject to constraints. 
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31 . (new): The method of Claim 26, wherein the wavelet transforms applying step 
b) comprises: 

folding the data matrix D into a (x+1)-dimensional multiway array D 
consisting of x spatial dimensions and 1 spectral dimension comprising p spectral 
channels, wherein x = 1, 2, or 3; 

applying an independent wavelet transform to each of the x spatial 
dimensions for each of the p spectral channels to provide a transformed multiway 

array D ; 

decimating the detail coefficients of the transformed multiway array D to 
obtain a multiway approximation array D a ; and 

unfolding the multiway approximation array D a to obtain the 
approximation matrix D a . 
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